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ABSTRACT 


Presently,  the  United  States  Navy  is  searching  for  an 
inqproved  method  to  predict  the  damage  to  a  ship  hull  or 
underwater  structure  that  results  from  an  underwater 
explosion.  One  method  of  predicting  this  damage  is  through 
the  use  of  nonlinear  finite  and  boundary  element  analysis. 
Underwater  Shock  Analysis  (USA)  combined  with  VEC/DYNA3D  is 
used  for  the  analysis  of  the  effect  of  explosive  shock  on 
numerical  models.  Initial  imperfections  are  introduced  in 
the  numerical  model  using  modal  imperfections.  The 
resulting  numerical  model  is  then  subjected  to  a  simulated 
underwater  shock  using  the  combined  USA/DYNA3D  code.  A 
sensitivity  analysis  is  performed  to  look  into  the  details 
on  the  damage  resulting  from  these  simulations. 
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1.  ZNTRODOCTION 


The  use  of  USA/DYNA3D  for  analyzing  structures 
subjected  to  underwater  shock  has  been  shown  to  be  an 
effective  tool  in  predicting  the  response  of  these 
structures  and  the  resulting  damage  (Chisum  1992).  In  some 
cases  the  actual  deformation  of  a  cylinder  subjected  to  an 
underwater  shock  is  not  accurately  predicted.  One 
explanation  for  this  inaccuracy  is  that  actual  cylinders 
have  many  imperfections  (e.g.,  out  of  roundness,  thin 
sections,  voids,  etc.)  and  the  finite  element  modeling  of 
these  cylinders  often  do  not  take  into  account  the  initial 
imperfections  that  are  present.  By  introducing 
inperfections  in  the  position  of  the  node  points  generated 
by  a  finite  element  mesh  generator  it  is  possible  to  more 
accurately  model  the  geometry  of  an  actual  cylinder.  With  a 
more  accurate  model  of  the  actual  cylinder  the  finite 
element  numerical  analysis  of  the  resulting  damage  due  to 
underwater  shock  may  be  more  accurate  as  well. 
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XI .  INITIAL  IMPBRFBCTIONS 


The  location  of  the  node  points  of  a  structure  modeled 
by  a  finite  element  mesh  generator  are  precisely  located  (to 
the  numerical  accuracy  of  the  modeling  code)  at  the 
positions  specified  in  the  inputs  to  the  mesh  generator 
algorithm.  Manufactured  structures  usually  have  many 
inperfections  that  are  not  accurately  modeled  by  many  finite 
element  mesh  generators.  Surv^s  have  been  performed  to 
measure  the  imperfections  that  naturally  exist  in 
cylindrical  shell  structures  (Arbocz  1982) .  These  data  have 
been  collected  into  data  banks,  and  Arbocz  suggests  that 
they  be  used  to  irrprove  design  criteria  for  *>uckling  of  thin 
shells.  The  inperfection  data  banks  show  that  inperfections 
have  characteristic  distributions  that  include  decreasing 
modal  anplitudes  with  increasing  mode  number.  Kirkpatrick 
(1989)  found  that  by  introducing  initial  modal  inperfections 
in  a  cylindrical  !=tructure  subjected  to  blast  shock  loading 
that  the  resulting  nxamerical  analysis  agreed  much  closer  to 
the  experimental  data  than  the  analysis  without  these 
initial  inperfections. 

Initial  inperfections  were  introduced  into  the  location 
of  the  node  points  of  the  finite  element  model  using  a 
siimmation  of  modal  inperfections  expressed  as  the  cosine 
series  shown  in  Eg.  (1) . 
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(1) 


N 

Ai?(0) »  J^AiC0s(«6+^) 

■>3 

where  ^  is  the  radial  iirqperfection,  6  is  the  angular 
position,  N  is  the  maximum  itKsdal  contribution,  A.  is  the 
modal  amplitude,  and  is  a  random  modal  phase  shift.  The 

assumption  that  the  modal  phase  shift  is  a  random  variable 
is  reasonable  for  many  shells.  An  enpirical  form  for  the 
modal  amplitude  is  expressed  as  shown  in  Eg.  (2). 


where  A«  is  the  modal  amplitude  of  the  nth  modal 
inperfection,  n  is  the  mode  number  and  X  and  r  are 
coefficients  used  to  fit  the  available  data  for  shells  of  a 
given  construction.  The  modal  anplitude  can  be  modeled  as  a 
constant  percentage  of  the  shell  thickness.  For  many  of  the 
studies  in  this  investigation  the  modal  amplitude  is 
constant  as  shown  in  Eg.  (3)  at  1%  of  the  shell  thickness 
for  all  modes. 

AnsQ.OlA  for  all  n . (3) 
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Kirkpatrick  found  that  the  modal  amplitudes  in  Eq.  (4)  and 
(5)  accurately  modeled  the  inperfections  of  his  test 


cylinder. 

i4ii«0,05/i  iiS6 . (4) 

n't!  . (5) 

n 


The  modal  imperfections  introduced  in  this  study  are 
assumed  to  be  a  function  only  of  the  circumferential 
position  of  the  node  and  not  a  function  of  the  axial 
position  of  the  node. 
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ZZZ.  TNO  DZMEMSZONAL  ZMFIMITB  CYLIHDBR  MODELS 


Fox  (1992)  developed  a  two  dimensional  infinite 
Q^iinder  model  to  validate  the  use  of  the  USA/DYNASD  code 
for  use  in  underwater  shoclc  analysis.  This  model  was 
further  developed  into  the  full  symmetry  two  dimensional 
infinite  cylinder  model  shown  in  Figure  1.  This  model  has 
40  elements  around  the  circumference  of  the  cylinder,  radius 
of  6  inches,  shell  thic)cness  of  0.06  inches  and  an  element 
length  of  0.006  inches.  The  same  restriction  on  the  length 
of  the  elements  found  for  the  model  Fox  developed 
(length/radius=0 . 001)  was  also  true  for  this  model.  Fox 
stated  that  longer  elements  resulted  in  'residual  three 
dimensional  effects'  that  caused  oscillations  in  the 
numerical  results.  Symmetric  boundary  conditions  were 
inqposed  on  both  ends  of  the  cylinder  to  model  an  infinite 
length  cylinder.  This  model  and  all  subsequent  models 
(unless  specifically  noted)  use  Belytsch)co-Tsay  (Belytschlco 
1984)  shell  elements  for  their  numerical  efficiency  and  to 
eliminate  the  problems  of  using  a  solid  element  to  model  a 
thin  walled  structure.  The  shell  material  is  modeled  as 
mild  steel,  an  elastic-plastic  material  with  a  Young's 
n^ulus  of  2.9x10*^  psi  ,  a  Poisson's  ratio  of  0.3,  a 
hardening  modulus  of  5100  psi  and  a  yield  stress  of  32,000 
psi. 


5 


6 


A.  F«rf«ct  Cyllndar  Modttl 

This  cylinder  was  subjected  to  a  underwater  shock  wave 
siodeled  as  a  plane  wave  with  an  anqplitude  of  500  psi  and  a 
duration  of  1  millisecond  resulting  from  an  explosive  charge 
located  along  the  x-axis.  The  early  damage  patteim  at  0.5 
milliseconds  is  shown  in  Figure  2.  The  front  face  of  the 
cylinder  is  flattened  by  the  impacting  plane  wave 
approaching  from  the  right.  However,  Figure  3  shows  that 
the  final  damage  pattern  at  5  milliseconds  is  very  different 
than  this  early  deformation  pattern.  All  plastic 
deformation  has  been  conqpleted  by  5  milliseconds.  The  side 
of  the  cylinder  facing  the  approaching  shock  wave  shows  a 
pronounced  protrusion.  This  pattern  was  unexpected  based  on 
experimental  results  which  have  usually  shown  inward 
deformations  of  a  cylindrical  shell  facing  the  explosive 
charge.  It  is  hypothesized  that  this  raised  section  is  due 
to  the  fact  that  the  modeled  cylinder  is  a  perfect  cylinder 
and  is  not  representative  of  a  true  manufactured  cylinder. 
The  manufactured  cylinder  will  have  many  imperfections  in 
the  shell  which  is  not  present  in  the  computer  models. 

In  order  to  study  the  effect  of  mesh  size  on  the 
response  of  this  cylinder  another  model  was  developed  with 
64  elex.  :its  around  the  circumference  of  the  cylinder. 
Figure  4  shows  that  the  response  of  this  cylinder  is  similar 
to  that  for  the  40  element  cylinder.  There  is  still  a 
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♦-Shock 

Wave 


0.5  ms  (deformations  scaled  up  lOOX) 

Figv.«  2.  Deformation  of  Two  Dimensional  Infinite  Cylinder 
Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for  1ms 


♦-Shock 

Wave 


5.0  ms  (deformations  scaled  up  lOX) 

Figure  3.  Deformation  of  Two  Dimensional  Infinite  Cylinder 
Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for  1ms 
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♦-Shock 

Wave 


5.0  ms  (deformations  scaled  up  lOX) 

Figure  4.  Deformation  of  Two  Dimensional  Infinite  Clinder 
Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for  ims 

Fine  Mesh 


distinct  protrusion  of  the  shell  toward  the  direction  of  the 
approaching  shock  wave.  All  subsequent  cylinders  will  use 
40  elements  around  the  circxunference  of  the  cylinder. 

B.  ZB^rfact  Cylinder  Models 

The  introduction  of  initial  inqper feet ions  significantly 
changes  both  the  shape  and  magnitude  of  the  resulting 
deformation.  Introduction  of  a  6th  mode  in^jer  feet  ion,  as 
shown  in  Figure  5  (imperfection  magnitude  scaled  by  a  factor 
of  100)  with  a  modal  airqplitude  of  5%  of  the  shell  thickness 
and  no  random  phase  shift,  results  in  the  damage  pattern 
shown  (deformations  scaled  lay  factor  of  10)  in  Figure  6.  A 
mode  6  inperfection  was  chosen  because  of  the  ease  in 
identifying  its  distinctive  pattern.  Use  of  another  model 
imperfection  would  only  have  changed  the  shape  of  the  final 
deformation  not  any  conclusions  resulting  from  this  study. 
Coitparison  of  the  final  damage  pattern  with  the  shape  of  the 
initial  imperfection  shows  that  the  damage  pattern  follows 
the  initial  inperfection  resulting  in  a  final  deunage  pattern 
looking  like  the  mode  6  initial  imperfection  magnified  due 
to  the  external  shock  wave  applied  to  the  cylinder.  The 
introduction  of  inperfections  caused  the  deformation  of  the 
cylindrical  shell  to  follow  the  initial  imperfection  pattern 
of  the  shell.  In  addition  to  changing  the  damage  pattern. 
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4-Shock 

Wave 


Initial  Mode  6  Iirqperfection  (imperfections  scaled  up  lOOX) 


Figure  5.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  1ms  Mode  6  Initial 

Inperfection  Ai^0.05h 


Figure  6.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  1ms  Mode  6  Initial 

Irrmerfection  A6  =  0.05h 


the  initial  inqperfection  also  causes  much  greater 
deflections  of  the  shell  for  the  same  shock  wave  intensity. 

Another  cylinder  was  modeled  with  inqper  feet  ions 
containing  the  first  10  modal  imperfections,  as  shown  in 
Figure  7  with  modal  amplitudes  for  each  mode  of  1%  of  the 
shell  thickness  and  with  random  phase  shifts.  The  resulting 
damage  pattern  is  shown  in  Figure  8.  Again,  the  general 
shape  of  the  final  deformation  pattern  is  very  close  to  the 
shape  of  the  initial  inperfection  pattern. 

The  introduction  of  a  different  initial  inperfection, 
shown  in  Figure  9  with  initial  imperfections  magnified  by  50 
times,  using  the  modal  eutplitude  of  equations  (4)  and  (5), 
results  in  the  final  'Reformation  pattern  shown  in  Figure  10. 
Again,  the  deformatic  of  the  shell  of  the  cylinder  follows 
the  initial  imperfection  pattern. 
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Initial  In^erfection  (in^er feet ions  scaled  up  200X) 

Figure  7.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  Ims  Initial 
Inperfection  First  10  Modes  i4*s0.01A  with  Random  Phase  Shift 


4-Shock 

Wave 


5.0  ms  (deformations  scaled  up  2X) 

Figure  8.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  1ms  Initial 
Inperfection  First  10  Modes  Am-O.Olh  with  Random  Phase  Shift 


3 


Initial  Imperfection  (imperfections  scaled  up  SOX) 


Figure  9.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  1ms  Initial  Inperfection 
First  10  modes  A«  =  (X05/i  for  Modes  2  through  6,  Am^lhln^  for 

Modes  7  through  10 


5.0  ms  (deformations  scaled  up  2X} 


4-Shock 

Wave 


Figure  10.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  1ms  Initial  Inperfection 
First  10  modes  Am^0.05h  for  Modes  2  through  6,  Am^lh/n^  for 

Modes  7  through  10 


ZV.  THIUn  DZMBNSIOHAL  RING  STIFPBMSD  IHFZMZTB  CTLINDER 


MODELS 


The  three-dimensional  ring  stiffened  infinite  cylinder 
model  shown  in  Figure  11  was  developed  from  the  two- 
dimensional  infinite  cylinder  model.  This  model  has  a  total 
of  400  elements  with  40  elements  in  the  circumferential 
direction  and  10  elements  along  the  length.  Again  symmetric 
boundary  conditions  were  in^osed  on  both  ends  of  the 
cylinder  to  model  an  infinite  cylinder.  This  three- 
dimensional  inodel  has  stiffeners  0.12  inches  thick  and  1 
inch  deep  located  on  12  inch  spacing.  The  shell  is  0.06 
inches  thick.  The  shell  and  stiffeners  are  modeled  as  mild 
steel  with  the  same  properties  as  listed  earlier 

A.  Perfect  Cylinder  Model 

This  cylinder  was  subjected  to  the  same  plane  shock 
wave  as  the  previous  models.  The  early  deformation  of  the 
cylinder  at  0.2  milliseconds  (Figure  12)  shows  the  outer 
shell  pinching  in  on  either  side  of  the  stiffener  located  in 
the  center  of  the  cylinder.  Figure  13  shows  the  flattening 
of  the  face  of  the  cylinder  facing  the  explosive  charge  at 
0.2  milliseconds.  At  this  time  in  the  deformation  history 
there  is  no  evidence  of  the  formation  of  a  local  raised  area 
on  the  front  of  the  cylinder  facing  the  charge.  The  final 
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Figure  11.  Three  Dimensional  Ring  Stiffened  Infinite 

Cylinder  Model 
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Front  View  0  0.2  ms  (deformations  scaled  up  SOX) 


Figure  12.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 

Pulse  of  500  psi  for  1ms 


<-Shock 

Wave 


End  View  0  0.2  ms  (deformations  scale  up  SOX) 

Figure  13.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 

Pulse  of  500  psi  for  1ms 


deformation  pattern  shown  in  Figure  14  shows  pinching  of  the 
shell  on  either  side  of  the  stiffener.  This  pinching  effect 
was  also  noted  by  Chisum  (1992)  and  appears  to  be  present  in 
many  finite  element  simulations  of  underwater  shock  which 
model  test  cylinders  as  perfect  cylinders.  Again,  the  final 
deformation  pattern  of  the  cylinder  side  facing  the  charge 
has  a  local  protrusion  of  the  shell  material  toward  the 
explosive  charge  as  shown  in  Figure  15.  Again,  it  is 
hypothesized  that  this  damage  pattern  is  the  result  of  the 
model  cylinder  being  a  perfect  cylinder  while  test  cylinders 
will  always  have  irrqper  feet  ions  in  the  shell  of  the  cylinder. 

B.  Zoqperfect  Cylinder  Models 

The  introduction  of  initial  imperfections  significantly 
changes  both  the  shape  and  magnitude  of  the  final 
deformations  seen  in  this  cylinder  model.  The  5%  mode  6 
in^erfection  shown  in  Figure  16  was  introduced  into  this 
model.  The  resulting  deformation  due  to  a  500  psi  1 
millisecond  plane  wave  pressure  pulse  is  shown  in  Figures  17 
through  20.  The  initial  imperfection  shape  is  clearly 
evident  in  the  final  damage  pattern.  However  the  shape  of 
the  cylinder  at  0.2  milliseconds  does  not  show  any  of  the 
initial  inperfection.  There  is  insufficient  time  elapsed 
for  the  deformations  to  grow  large  enough  to  show  this 
initial  inperfection.  However  the  pinching  of  the  shell  on 
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Front  View  0  5 . 0  ms 


Figure  14.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 

Pulse  of  500  psi  for  1ms 


<-Shock 

Wave 


End  View 

Figure  15.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 

Pulse  of  500  psi  for  1ms 


Mode  6  Initial  Imperfection  A«s0.05n  (imperfections  scaled 

up  lOOX) 

Figure  16.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 

Pulse  of  500  psi  for  1ms 
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Front  View  @  0.2  ms  (deformations  scaled  up  SOX) 


Figure  17.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  with  a  Mode  6  Initial  Imperfection 
At^O.OSh  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 

1ms 


<-Shock 

Wave 


End  View  d  0.2  ms  (deformations  scaled  up  SOX) 

Figure  18.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  with  a  Mode  6  Initial  In^erfection 
i46*0.05A  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 

1ms 


Front  View  @  5.0  ms 


Figure  19.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  with  a  Mode  6  Initial  Imperfection 
i4«  =  0.05/i  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 

1ms 


<-Shock 

Wave 


End  View  9  5.0  ms 

Figure  20.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  with  a  Mode  6  Initial  Inperfection 
AisiO.OSh  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 

1ms 


either  side  of  the  stiffener  has  already  begxin  at  0.2 
xtiilliseconds.  The  final  deformation  still  shows  this 
pinching  effect  however  the  magnitude  of  this  type  of 
deformation  is  reduced  by  the  initial  inperfection. 

With  the  introduction  of  an  initial  inperfection  of  the 
first  10  modes  with  modal  anplitudes  of  1%  of  the  shell 
thickness  and  with  random  phase  shifts  as  shown  in  Figure  21 
the  resulting  damage  pattern  again  shows  that  the 
deformation  of  the  shell  will  preferentially  follow  the 
initial  inperfections.  The  resulting  damage  pattern  is  very 
different  from  the  perfect  case  (Figures  22  and  23).  Most 
noteworthy  is  the  elimination  of  the  pinch  in  the  outer 
shell  on  either  side  of  the  stiffener.  In  addition  the 
magnitude  of  the  deformation  in  the  outer  shell  is  much 
greater  than  for  the  perfect  cylinder  case.  The 
introduction  of  these  inperfections  results  in  sites  where 
the  cylinder  preferentially  deforms  during  the  shock 
pressure.  Thus  when  the  cylinder  deforms  the  deformation 
follows  the  initial  inperfections  resulting  in  the  final 
shape  of  the  cylinder  looking  like  the  initial  inperfection 
shape. 
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Initial  Ixtqperfection  (imperfections  scaled  up  200X) 

Figure  21.  Three  Dimensional  Ring  Stiffened  Infinite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1ms  Initial  Inperfection  First  10  Modes  i4«  =  0.01A  with  Random 

Phase  Shift 


Front  View  0  5.0  ms 


Figure  22.  Three  Dimensional  Ring  Stiffened  Infinite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1ms  Initial  Imperfection  First  10  Modes  i4«  =  0.01A  with  Random 

Phase  Shift 


<-Shock 

Wave 


End  View  9  5 . .  0  ms 

Figure  23.  Three  Dimensional  Ring  Stiffened  Infinite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1ms  Initial  Imperfection  First  10  Modes  i4iis0.01A  with  Random 

Phase  Shift 


V.  TBXII  DZMBMSIONAL  RIMO  STZFFIMSD  FZNZTI  LENGTH  CYLZNDBR 


From  the  three-dimensional  infinite  cylinder  a  three 
dimensional  finite  length  cylinder  was  developed  as  shown  in 
Figure  24.  This  cylinder  is  three  feet  long,  1  foot  in 
diameter  with  two  stiffeners  evenly  spaced  12  inches  apart. 
The  shell  of  the  cylinder  is  mild  steel  0.06  inches  thick. 
The  stiffeners  are  mild  steel  0.12  inches  thick  and  1  inch 
deep.  The  endplate  is  HY-100  steel  0.25  inches  thick.  The 
HY-100  steel  is  modeled  as  an  elastic-plast  material  with 
a  Young's  modulus  of  2.9x10*^  psi,  a  Poisson's  ratio  of  0.3, 
a  hardening  modulus  of  5100  psi  and  a  yield  stress  of 
108,000  psi.  This  cylinder  is  modeled  as  a  half  cylinder 
with  a  plane  of  symmetry  perpendicular  to  the  axis  of 
rotation.  The  model  has  40  elements  in  the  circxunferential 
direction  and  15  elements  in  the  axial  direction  for  a  total 
of  600  elements  and  921  nodes.  The  use  of  a  half  symmetry 
model  with  a  symmetric  boundary  condition  results  in  a 
smaller  number  of  elements  and  greater  conputational 
efficiency  for  the  finite  element  analysis. 

A.  Perfect  Cylinder  Models 

This  cylinder  was  subjected  to  the  same  planar  shock 
wave  as  in  the  previous  models.  The  resulting  damage 
pattern  for  a  perfect  cylinder  can  be  seen  in  Figures  25 
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through  32.  Again  the  side  of  the  cylinder  facing  the 
planar  shock  wave  shown  a  local  raised  area  that  was  also 
seen  in  the  previous  models.  In  addition  there  is  a  pinch 
in  the  shell  near  the  endplates  and  on  either  side  of  the 
stiffeners.  This  pinch  was  also  noted  by  Chisum  (1992) 
during  his  analysis  of  DNA  models  used  for  underwater  shock 
testing.  Most  of  the  deformation  of  the  cylinder  shell 
occurs  between  1.0  and  2.0  milliseconds.  During  the  time 
from  1.0  to  2.0  milliseconds  this  kinetic  energy  is 
transferred  to  strain  energy  causing  deformation  of  the 
shell.  After  2.0  milliseconds  the  most  of  the  plastic 
deformation  has  occurred  and  the  cylinder  shell  then  deforms 
only  through  elastic  vibration.  Figure  33  is  a  plot  of  the 
kinetic  energy  of  the  cylinder  shell  as  a  function  of  time. 
The  shell  has  its  maximum  kinetic  energy  at  1  ms  and  then 
this  energy  rapidly  decreases  as  it  is  dissipated  into  the 
surrounding  medium  and  into  the  shell  material  as  strain 
energy.  Figure  34  shows  the  strain  history  of  the  element 
on  the  front  of  the  cylinder  exactly  in  the  middle  of  the 
stiffeners.  Again,  this  plot  shows  that  most  of  the  hoop 
strain  occurs  between  1  and  2  milliseconds. 

B.  Ziqperfect  Cylinder  Models 

The  addition  of  initial  inperfections  greatly  changes 
the  resulting  deformation  pattern  of  the  shell  of  the 
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Front  View  d  0.2  ms  (deformations  scaled  up  SOX) 

Figure  25.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 


Shock-> 

Wave 


End  View  9  0.2  ms  (deformations  scaled  up  50X) 

Figure  26.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 
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Front  View  d  1.0  ms  (deformations  scaled  up  2X) 

Figure  27.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 


End  View  d  1.0  ms  (deformations  scaled  up  2X) 

Figure  28.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 
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Front  View  d  2.0  ms  (deformations  scaled  up  2X) 

Figure  29.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 


Shock-> 

Wave 


End  View  @  2.0  ms  (deformations  scaled  up  2X) 

Figure  30.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 
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Front  View  9  5.0  ms  (deformations  scaled  up  2X) 


Figure  31.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 


Shock-^ 

Wave 


End  View  9  5.0  ms  (deformations  scaled  up  2X) 

Figure  32.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 


cylinder.  The  addition  of  a  node  6  iirperfection  (shown  in 
Figure  35)  with  a  modal  an?)lituc3  of  5%  of  the  shell 
thickness  and  with  no  random  phase  shift  results  in  the 
final  damage  pattern  shown  in  Figure  36  and  37.  The  end 
view  clearly  shows  the  strong  effect  of  the  mode  6  initial 
imperfection  causing  the  final  deformation  to  follow  this 
initial  imperfection.  The  front  view  of  the  cylinder  shows 
that  the  pinching  of  the  shell  near  the  endplates  and 

stiffeners  is  very  much  reduced  due  to  the  introduction  of 
this  inperfection. 

An  investigation  into  the  effect  of  the  modal  amplitude 
on  the  final  deformation  pattern  showed  that  very  small 
modal  anplitudes  have  a  significant  effect  o  the  response 
of  the  cylinder.  An  initial  inqperfection  consisting  of  the 
first  10  mode  shapes  with  varying  modal  amplitudes  and 

random  phase  shifts  is  shown  in  Figure  38.  The  resulting 
damage  patterns  at  various  modal  anplitudes  are  shown  in 
Figures  39  through  50.  The  shell  deformation  again  follows 
the  initial  imperfections  present  in  the  shell  structure. 
As  the  modal  anplitudes  increase  the  deformation  of  the 
shell  also  increases.  At  a  modal  amplitude  of  0.01%  of  the 
shell  thickness  the  front  of  the  cylinder  is  distinctly 
raised  toward  the  explosive  charge  as  in  the  perfect 
cylinder.  But,  it  can  also  be  seen  from  the  end  view  that 

the  deformation  pattern  even  at  this  magnitude  of 
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Initial  Ixrqper  feet  ion  (in^erfections  scaled  up  lOOX) 


Figure  35.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
Ims  Mode  6  Initial  In5>erfection  Ae^^O.OSh 


Front  Viev»  9  5.0  ms 

Figure  36.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1ms  Mode  6  Initial  litter  feet  ion  A6=s0.05/i 


Shock-4 

Wave 


End  View  9  5.0  ms 

Figure  37.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1ms  Mode  6  Initial  iR^erf action  A6s:0.0SA 
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Initial  Iirperf action  (imperfections  scaled  up  200X) 

Figure  38.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Sxibjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1ms  First  10  Modes  w/Random  Phase  Shift 


Front  View  A«s0.0001A 


Figure  39.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 


End  View  i4iis0.0001A 


Figure  40.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Sxibjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 


Front  View  i4>s0.000SA 


Figure  41.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  Ixns 


Shock-> 

Wave 


End  View  i4-  =  0.0005A 


Figure  42.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1ms 


Front  View  Ams0.00\h 


Fi^re  43.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
Ims  Deformations  at  5  milliseconds 


Shock“-> 

Wave 


End  View  A.^O-OOIA 

Figure  44.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1ms  Deformations  at  5  milliseconds 


Front  View  i4«*0.005A 


Figure  45.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  S^are  Pressure  Pulse  of  500  psi  for 
1ms  Deformations  at  5  milliseconds 


Shock-4 

Wave 


End  View  Am^O.OOSh 

Figure  46.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1ms  Deformations  at  5  milliseconds 


Front  View  A«sb0.01A 


Fi^re  47.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
Ims  Deformations  at  5  milliseconds 


Shock-» 

Wave 


End  View  i4as0.01A 

Figure  48.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  S^are  Pressure  Pulse  of  500  psi  for 
1ms  Deformations  at  5  milliseconds 


Front  View  AukO.OSA 

Figure  49.  Iliree  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  S^are  Pressure  Pulse  of  500  psi  for 
Ims  Oefoimations  at  5  milliseconds 


Shock-> 

Wave 


End  View  Am^0.05h 

Figure  50.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  S<^are  Pressure  Pulse  of  500  psi  for 
1ms  Deformations  at  5  milliseconds 
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inperfection  has  changed  from  the  pattern  for  the  perfect 
cylinder.  At  a  model  anqplitude  of  0.05%  this  alteration  of 
the  pattern  is  quite  obvious.  The  pinching  of  the  shell  of 
the  cylinder  persists  until  higher  modal  anplitudes  are 
reached.  Finally  at  a  modal  amplitude  of  0.5  %  the  pinching 
of  the  shell  is  eliminated.  It  should  be  noted  that  an 
amplitude  of  0.5%  of  the  shell  thickness  represents  a  modal 
amplitude  of  0.0003  inches  for  this  model.  This  magnitude 
of  imperfections  would  almost  certainly  be  present  in  a  test 
cylinder. 

The  development  of  a  raised  section  on  the  cylinder 
facing  the  shock  wave  was  unexpected.  With  the  previous 
perfect  cylinder  the  shock  wave  impacted  the  cylinder  along 
the  line  of  nodes  at  the  front  of  the  cylinder.  The 
cylinder  was  then  rotated  so  that  the  shock  wave  would 
impact  the  cylinder  along  a  line  of  elements.  The  resulting 
deformation  pattern  is  shown  in  Figure  51.  There  is  still  a 
local  raised  section  on  the  cylinder  facing  the  shock  wave 
direction.  Thus,  it  is  felt  that  the  shot  geometry  of  the 
finite  element  mesh  is  not  a  contributing  factor  to  the 
development  of  this  raised  section. 

In  addition,  a  very  specific  inperfection  was 
introduced  in  the  five  nodes  facing  the  explosive  charge 
centrally  located  in  the  midbay  of  the  cylinder.  This 
imperfection  was  an  inward  imperfection  of  the  node 


44 


Shock-» 

Wave 


Figure  51.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
Ims  Cylinder  Rotated  4.5  Degrees  End  View  Final  Deformation 

at  5  milliseconds 
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locations  with  an  anplitude  of  1%  of  the  shell  thickness 
(0.006  inches).  The  damage  pattern  resulting  from  a  500  psi 
1  millisecond  shock  wave  is  shown  in  Figure  52. 

The  use  of  Hughes-Liu  shell  element  (Hughes  1981)  vice 
the  Belytschko-Tsay  shell  element  was  investigated.  The 
Hughes-Liu  shell  element  does  not  use  many  of  the 
sisplifying  assunptions  used  in  formulating  the  Belytschko- 
Tsay  shell  element.  As  a  result «  the  use  of  the  Hughes-Liu 
shell  element  requires  longer  computational  times  that  may 
not  be  acceptable  in  certain  models.  Figures  53  and  54  show 
the  resulting  deformation  to  a  model  using  Hughes-Liu  shell 
elements.  This  model  is  identical  to  the  previous  models 
with  the  exception  of  the  use  of  the  Hughes-Liu  shell 
element.  This  model  has  initial  imperfections  with  the 
first  10  modes  and  modal  amplitudes  of  1%  of  the  shell 
thickness  and  random  phase  shifts.  Comparison  of  Figures  53 
and  54  with  Figures  47  and  48  shows  that  the  Hughes-Liu  and 
the  Belytschko-Tsay  shell  el«nents  have  the  same  final 
deformations.  In  this  case  the  use  of  the  niunerically  more 
complicated  Hughes-Liu  shell  element  is  not  warranted. 

C.  Exponential  Decay  Shock  Wave 

The  use  of  an  explosive  shock  wave  modeled  as  a  square 
wave  is  representative  of  a  very  large  explosion  at  a  great 
standoff  distance  from  the  cylinder.  In  order  to  model  a 
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Figure  52.  Three  Dimensional  Finite  cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  litis  End  View  Initial 
Imperfection  Front  Midbay  Nodes  1%  Offset  Radially  Inward 
Final  Deformation  at  5  milliseconds 
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Front  View  9  5.0  ms 

Figure  53.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Hughes-Liu  Shell  Elements  Subjected  to  a  Square 
Pressure  Pulse  of  500  psi  for  1  ms  Initial  Imperfection 
First  10  Modes  i4.s0.01A  with  Random  Phase  Shift  Final 
Deformation  at  5  milliseconds 


Shock-» 

Wave 


End  View  9  5.0  ms 

Figure  54.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Hughes-Liu  Shell  Elements  Subjected  to  a  Square 
Pressure  Pulse  of  500  psi  for  1  ms  Initial  Imperfection 
First  10  Modes  i4iia0.01A  with  Random  Phase  Shift  Final 
Deformation  at  5  milliseconds 
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smaller  explosion  close  to  the  cylinder  a  pressure  profile 
corresponding  to  that  produced  by  40  pounds  of  PETN  at  a 
standoff  distance  of  30  feet  was  used.  This  explosion 
produces  a  peak  pressure  of  1839  psi  and  an  exponentially 
decaying  pressure  history  with  a  time  constant  of  0.377 
milliseconds.  The  resulting  deunage  to  the  cylinder  from 
this  shock  wave  is  shown  in  Figure  55  through  57.  From  the 
offset  view  it  can  be  seen  that  the  response  of  the  cylinder 
is  different  from  that  to  the  previous  plane  wave.  The 
small  standoff  distance  causes  the  shock  wave  to  be  a 
spherical  wave  at  the  cylinder.  The  midbay  still  has  a 
protrusion  of  the  shell  material  toward  the  explosive  charge 
however  the  endbay  shell  material  has  a  depression  of  the 
shell  facing  the  charge.  Also  evident  from  the  front  view 
is  the  severe  pinching  of  the  shell  near  the  endplate  and  on 
either  side  of  the  stiffener.  This  pinching  was  very  severe 
near  the  endplate. 

The  introduction  of  a  5%  mode  6  iitperfection  to  the 
model  (Figure  58)  results  in  the  deformations  shown  in 
Figure  59  and  60.  Again,  it  is  evident  that  the  damage 
pattern  of  the  cylinder  due  to  the  explosive  shock  clearly 
followed  the  initial  inperfection.  The  pinching  of  the 
shell  material  near  the  endplate  and  stiffeners  was  reduced 
but  not  eliminated. 
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Front  View 


Figure  55.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Final  Deformation  at  5  milliseconds 


Offset  View 

Figure  56.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Final  Deformation  at  5  milliseconds 


Shock-4 

Wave 


End  View 

Figure  57.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Final  Deformation  at  5  milliseconds 
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Initial  Imperfection  (imperfections  scaled  up  lOOX) 

Figure  58.  Three  Dimensional  Ring  Stiffened  Finite  Length 
cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Initial  Mode  6  Imperfection  Ai-O.OSh 


Front  View  9  5.0  ns 

Figure  59.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Initial  Mode  6  Imperfection  i46s0.05A 


Shock-» 

Wave 


End  View  @  5.0  ms 

Figure  60.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Initial  Mode  6  Inperfection  A6-0.0Sh 
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Introducing  an  inperf action  of  the  first  10  mode  shapes 
with  a  modal  amplitude  of  1%  of  the  shell  thickness  with 
random  i^se  shifts  (Figure  61)  results  in  the  final  damage 
pattern  shown  in  Figure  62  and  63.  The  deformation  of  the 
shell  of  the  cylinder  was  again  influenced  by  the  initial 
imperfection  pattern.  The  local  pinching  of  the  shell 
material  on  either  side  of  the  stiffener  was  eliminated  and 
the  pinching  near  the  endplate  was  reduced  in  magnitude. 


Initial  Imperfection  (imperfections  scaled  up  200X) 

Figure  61.  ^ree  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Initial  Imperfection  First  10  inodes  >1-.  =  0.01A 

with  Random  Phase  Shift 
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Front  View  9  5.0  xns 


Figure  62.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Initial  Imperfection  First  10  modes  >4ii  =  0.01/i 

with  Random  Phase  Shift 


End  View  9  5.0  ms 

Figure  63.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff  Initial  Imperfection  First  10  modes  A»  =  0.0lh 

with  Random  Phase  Shift 
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VZ.  SUmUJlT  AMD  C0MCLU8I0MS 


The  response  of  test  cylinders  to  underwater  explosive 
shock  is  a  explicated  function  of  many  factors.  One 
important  factor  that  is  often  overlooked  in  the  modeling  of 
these  underwater  shock  phenxiena  is  the  initial 
imperfections  that  are  present  in  the  test  cylinders.  The 
response  of  model  cylinders  subjected  to  simulated 
underwater  explosive  shock  has  been  shown  to  be  very 
dependent  to  the  initial  inperfections  introduced  into  the 
cylinder  models.  The  introduction  of  these  inperfections 
not  only  caused  the  shape  of  the  shell  of  the  cylinder  (as 
viewed  from  the  end)  to  follow  the  shape  of  the  initial 
imperfection  but  it  has  also  changed  the  response  of  the 
shell  near  the  endplates  and  stiffeners  of  the  model 
cylinder.  The  pinching  of  the  shell  near  these  stiffeners 
and  endplates  was  greatly  reduced  or  eliminated  by  the 
introduction  of  initial  inperfections  in  the  model 
cylinders.  This  resulting  response  is  closer  to  the 
response  observed  in  test  cylinders  subjected  to  actual 
underwater  shock  loading.  In  addition,  the  initial 
imperfections  introduced  increased  the  magnitude  of  the 
shell  deformations  conpared  to  the  deformations  for 
cylinders  modeled  as  perfect  cylinders.  If  the  initial 
imperfections  of  a  test  cylinder  are  known  and  introduced 


56 


into  a  modal  cylinder  used  in  finite  element  analysis,  the 
results  of  the  finite  element  analysis  may  simulate  the 
actual  response  of  the  test  cylinder  more  closely. 
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APPIMDZX  A 


INGRID  INPUT  FILE  FOR  TWO  DIMENSIONAL  INFINITE  CYLINDER 

MODEL 


INFINITE  CYLINDER  MODEL 

dn3d  vec  term  5.0e-3  plti  l,0e-5  prti  1000. 

mat  1  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e-«-4  sigy  3.2e*»-4  shell  quad  5  thick  0.060  endmat 

led  1  2  0.0  0.0  1.0  0.0 
led  2  2  0.0  2.5e-8  1.0  2.5e-8 


plan  2 

0  0  -.003  0  0-1  0.001  syiren 

0  0  .003  0  0  1  0.001  symm 

start 

-1  6  -11  ; 

-1  6  -11  ; 

1  2  ; 

-1.  0.  1. 

-1.  0.  1. 

-.003  .003 

all033036.0 
pri  -1-3  ;  -1-3  ;  ;  1  -1.0  0.  0.  0. 
mate  1 
end 
end 
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APPEMDIX  B 


INGRID  INPUT  FILE  FOR  THREE  DIMENSIONAL  STIFFENED  INFINITE 

CYLINDER  MODEL 


STIFFENED  INFINITE  CYLINDER  MODEL 

dn3d  vec  term  5.0e-3  plti  l.Oe-5  prti  1000. 

c 

c  material  definitions 

c 

c  shell  material  (mild  steel) 

c 

mat  1  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e-»-4  sigy  3.2eH-4  shell  quad  5  thiclc  0.060  endmat 
c 

c  stiffener  material  (mild  steel) 

c 

mat  2  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e4-4  sigy  3.2e'»-4  shell  quad  5  thiclc  0.12  endmat 

c 

c  load  curve  definitions 
c 

led  1  2  0.0  0.0  1.0  0.0 
led  2  2  0.0  l.Oe-6  1.0  l.Oe-6 

c 

c  symmetry  planes 

c 

plan  2 

0  0  -6.  0  0-1  0.001  symm 

0  0  6.  0  0  1  0.001  symm 

c 

c  construct  shell 

c 

start 

-1  11  -21  ; 

-1  11  -21  ; 

1  6  11  ; 

-1.  0.  1. 

-1.  0.  1. 

—6 .  0 .  6 . 
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all033036.0 
pri  -1-3  ;  -1-3  ;  ;  1  -1.0  0.  0.  0. 
mate  1 
end 
c 

c  construct  stiffener 
c 

start 

1  2  12  22  23  ; 

1  2  12  22  23  ; 

-1  ; 

-1-1011 

-1-1011 

0 

di  1  2  0  4  5  ;  1  2  0  4  5  ;  ; 

alio  550  3  6.0 

a220  440  3  5.0 

d  2  2  0  4  4  0 

mate  2 

end 

end 
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APPENDIX  C 


INGRID  INPUT  FILE  FOR  THREE  DIMENSIONAL  STIFFENED  FINITE 

CYLINDER  MODEL 


STIFFENED  FINITE  CYLINDER 

dn3d  vec  term  5.0e-3  plti  l.Oe-5  prti  1000. 

c 

c  material  definitions 

c 

c  shell  material  (mild  steel) 

c 

mat  1  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e-t>4  sigy  3.2e+4  shell  quad  5  thic)c  0.060  endmat 
c 

c  stiffener  material  (mild  steel) 

c 

mat  2  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e+4  sigy  3.2e+4  shell  quad  5  thick  0.12  endmat 

c 

c  endplate  material  (HY-lOO  steel) 

c 

mat  3  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.02e'»'4  sigy  1.08e+5  shell  quad  5  thick  0.25 
endmat 
c 

c  load  curve  definitions 
c 

led  1  2  0.0  0.0  1.0  0.0 
led  2  2  0.0  l.Oe-6  1.0  l.Oe-6 

c 

c  symmetry  planes 

c 

plan  1 

0  0  0.0  0  0  -1  0.001  symm 

c 

c  construct  shell 

c 
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start 

-1  6  -11  ; 

-1  6  -11  ; 

1  16  ; 

-1.  0.  1. 

-1.  0.  1. 

0.0  18.0 
all033036.0 
PX’i  -1-3  ;  -1-3  ;  ;  1  -1.0  0.  0.  0. 
mate  1 
end 


c 

c  construct  stiffener 
c 


start 


1271213; 

1  2  7  12  13  ; 

-1  ; 

-1-1011 

-1-1011 

6.0 

di  1  2  0  4  5  ;  1  2  0  4  5  ;  ; 

alio  550  3  6.0 

a220  440  3  5.0 

d  2  2  0  4  4  0 

mate  2 

end 


c 

c  construct  endplate 

c 

c  surface  definitions 

c 

sd  1  cyli  000  001  6.0 

sd  2  cyli  000  001  [6. 0*3/5] 

start 

1  5  10  15  19  ; 

1  5  10  15  19  ; 

-1  ; 

-1-1011 

-1-1011 

18.0 

prill  551  1  -1.0  0.0  0.0  0.0 

di  1  2  0  4  5  ;  1  2  0  4  5  ;  ; 

sfvi  -2  -4  ;  -2  -4  ;  ;  sd  2 

Sfi  -1  -5  ;  -1  -5  ;  ;  sd  1 

mate  3 
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APPENDIX  D 


FORTRAN  PROGRAMS  FOR  MODIFYING  INGRIDO  FILE  FOR  MODAL 

IMPERFECTIONS 


program  imp 

*  This  program  written  by  Donald  T.  Hooker  II  on  3/15/93  * 

*  Revised  5/12/93  * 


*  This  program  modifies  the  nodes  of  a  file  called  * 

*  "ingrido.raw" ,  which  is  an  edited  version  of  ingrido  in  * 

*  which  the  first  line  is  the  beginning  of  the  node  list.  * 

*  This  is  one  of  a  series  of  2  programs,  which  need  to  be  * 

*  run  in  order,  i.e.  "inp*  then  “convert*.  These  programs* 

*  modify  the  node  positions  for  modal  imperfections  WITH  * 

*  OR  WITHOUT  random  anplitudes  and/or  phase  shifts.  The  * 

*  final  output  file  name  is  inp.dat.  * 


* 

* 

* 

* 

* 

* 

* 

* 

* 


* 

* 

* 

* 


* 

* 

* 

** 


VARIABLE  DECLARATION 

a  -  angle  from  reference  axis  to  radial  position  of 
node 

al  thru  alO  -  weighting  factor  for  each  modal 
imperfection 
bl  -  ingrido  constants 
b2  -  ingrido  constants 

dr  -  change  in  radial  position  due  to  modal 
inperfections 

drl  thru  drlO  -  radial  position  change  for  each  modal 
imperfection 
h  -  thickness  of  shell 
n  -  node  number 
m  -  number  of  nodes 
pi  -  3.14159265359 

psl  thru  pslO  -  phase  shift  for  each  modal 
inperfection 

r  “  radial  position  of  node 

randanp  -  if  random  anplitude  is  desired  set  randanp=l 
randphs  -  if  random  phase  is  desired  set  randphs=l 
X  -  X  coordinate  of  node 
y  -  y  coordinate  of  node 
2-2  coordinate  of  node 

***♦*****♦****♦****♦**♦*'**♦♦♦****♦****♦*******♦♦****★★*★★ 


double  precision  x,y,z,a,r,dr,h,pi 
double  precision  al,a2,a3,a4,a5,a6,a7,a8,a9,al0 
double  precision  drl,dr2,dr3,dr4,dr5,dr6,dr7,dr8,dr9, 
&  drlO 

double  precision  psl,ps2,ps3,ps4,ps5,ps6,ps7,ps8,ps9, 
&  pslO 
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integer  randanp, randphs 


OPEN  INPUT  AND  OUTPUT  FILES 


open (14, files 'ingrido. raw* ) 
open (15, files ' imp. unconverted* ) 


IF  RANDOM  MODE  AMPLITUDE  IS  DESIRED  SET  RANDAMP  EQUAL* 
TO  1  * 

IF  RANDOM  MODE  PHASE  SHIFT  IS  DESIRED  SET  RANDPHS  * 
EQUAL  TO  1  * 


randamp  s  o 
randphs  s  i 

A**************'******************'***********'*********'******* 

*  INPUT  NUMBER  OF  NODES  TO  BE  MODIFIED  AND  WRITE  NUMBER* 

*  TO  OUTPUT  FILE  IMP . UNCONVERTED  * 


ms  921 
write(15, 98)m 
98  format (lx, i5) 

A*********#*'**'******************#***'**********************'** 

*  INPUT  THICKNESS  OF  SHELL  * 

*♦♦*♦****★****■»★*♦**♦★*★■*•♦**★■***★**♦**♦*****♦•****♦***♦♦*** 


h  s  0.06 

*  SET  THE  VALUE  OF  PI  * 


pi  =  3.14159265359 

***«****************«*««*4r«*****««******«********<lr***«r#**«^** 

*  INPUT  THE  WEIGHTING  COEFFICIENTS  OF  THE  FIRST  TEN  * 

*  MODE  SHAPES  * 

*****'*«'*******«'*'*«*«**'*******tlr*«***<**************#**«'***«.,^*# 


al  s  0.01 
a2  s  0.01 
a3  s  0.01 
a4  s  0.01 
a5  s  0.01 
a6  s  0.01 
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a7  «  0.49 
a8  «  0.64 
a9  >  0.81 
alO  «  1.0 


INPUT  PHASE  SHIFT  FOR  EACH  MODE  SHAPE 


psl  «  0.0 
ps2  s  0.0 
ps3  a  0.0 
ps4  s  0.0 
ps5  s  0.0 
ps6  s  0.0 
ps7  a  0.0 
ps8  a  0.0 
ps9  a  0.0 
psio  a  0.0 


*  INPUT  SEED  NUMBER  FOR  RANDOM  NUMBER  GENERATOR  * 

****************»*******************«****'*****'*****tlr******** 


seed  a  1.0 

*  CALL  SUBROUTINE  SRAND  TO  GENERATE  SEED  NUMBER  FOR  * 

*  RANDOM  NUMBER  GENERATOR  SUBROUTINE  * 


call  srand(seed) 


■trit********************************************************** 

*  CALL  SUBROUTINE  RAND  TO  GENERATE  RANDOM  PHASE  SHIFT  * 

*  FOR  EACH  MODE  SHAPE  * 

*  ♦********•♦*****♦■*♦***********♦♦***♦*♦*♦♦♦**•*♦♦*★•**★****♦* 


if  (randphs 
psl  a  2.  ♦ 
ps2  a  2.  * 
ps3  a  2.  * 
ps4  a  2.  * 
ps5  a  2.  * 
ps6  a  2.  * 
ps7  a  2.  * 
ps8  a  2.  * 
ps9  a  2 .  * 
PSlO  a  2, 
endif 


. eq . 1 )  then 
pi  *  randO 
pi  ♦  randO 
pi  *  rand() 
pi  *  rand() 
pi  *  rand() 
pi  *  randO 
pi  *  randO 
pi  *  randO 
pi  *  randO 
♦pi  *  randO 
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CALL  SUBROUTINE  RAND  TO  GENERATE  RANDOM  WEIGHTING  OF  * 
MODE  SHAPES  * 


if  (randeunp  .eq.  1)  then 
al  s  al  *  randO 
a2  s  a2  *  rand() 
a3  s  a3  *  rand(i 
a4  s  a4  *  randO 
a5  »  a5  *  rand() 
a6  s  a6  *  rand() 
a7  a  a7  *  rand ( ) 
a8  a  a8  *  rand ( ) 
a9  =  a9  *  randh 
alO  a  alO  *  randO 
endif 


A*********************'**'*****'****************************'*** 

*  DO  LOOP  TO  READ  INPUT  DATA  * 


do  100  isl,jn 

read(14,*)  n,bl,x,y, 2,b2 

A******'*'*'*'**'*'******'*'*'**'****'*'*'**'***************************'** 

*  CALCULATE  RADIUS  OF  NODE  POSITION  * 

*  (Look  carefully  at  the  axis  of  rotation  of  the  * 

*  cylinder  to  determine  the  proper  cartesian  * 

*  coordinates  used  to  calculate  the  radius  of  the  * 

*  cylinder.  In  this  case  the  axis  of  rotation  is  the  z* 

*  axis.  * 


ra (x*x+y*y ) **0.5 

*  CALCULATE  THE  ANGLE  FROM  THE  REFERENCE  AXIS  TO  THE  * 

*  NODE  POINTS  (Look  carefully  at  the  zero  angle  * 

*  direction.  In  this  case  the  x  eocis  points  in  the  * 

*  direction  of  the  charge.)  * 

if  ( (dabs (y) .It.  0.001))  then 
if  (x.gt.0.0)  then 
aaO.  0 
else 
a=pi 
endif 
else 

a=datan  (y/x) 
if  (x. It. 0.0)  then 


» 
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ava-^pi 

endif 

andif 


MODIFY  THE  NODE  POSITIONS  USING  THE  FIRST  TEN  MODE 
SHAPES 


& 


drl 

S 

al 

if 

h 

it 

dcos  (a 

+  psl) 

dr2 

s 

a2 

* 

h 

* 

dcos (2. 

*  a  +  ps2) 

dr3 

m 

a3 

h 

* 

dcos (3. 

*  a  +  ps3) 

dr4 

s 

a4 

* 

h 

* 

dcos (4. 

*  a  +  ps4) 

drS 

s 

a5 

h 

* 

dcos ( 5 . 

*  a  +  ps5) 

dr6 

s 

a€ 

* 

h 

* 

dcos ( 6 . 

*  a  +  ps6) 

dr7 

s 

a7 

it 

h 

/ 

7.  /  7. 

*  dcos (7.  * 

a 

+ 

ps7) 

dr8 

= 

a8 

it 

h 

/ 

8.  /  8. 

*  dcos ( 8 .  * 

a 

+ 

ps8) 

dr9 

s 

a9 

it 

h 

/ 

9.  /  9. 

*  dcos ( 9 .  * 

a 

+ 

ps9) 

drlO 

8 

:  alO 

* 

h 

/  10.  / 

10.  *  dcos (10 

* 

a  +  pslO) 

dr  = 

drl 

+ 

dr2 

+  dr3  + 

dr4  +  dr5  + 

dr6 

+  dr7  +  dr8  + 

dr 9  +  dr 10 


*  CALCULATE  THE  CARTESIAN  COORDINATES  OF  THE  NODES  * 

A*'***'******'***"**'*'**'******'****'************'************'******* 


xsx--dr*dcos  (a) 
y»y-dr*dsin(a) 

************************************************************* 

*  WRITE  DATA  TO  TEMPORARY  DATA  FILE  IMP .  UNCONVERTED  * 

***********'***********'************«****«*********«««***«*«** 

write(15, 998)  n,bl,x,y, z,b2 
998  fonnat(i8,f5.0,3e20.13,f5,0) 

100  continue 

*'*‘*'*‘******«*«r'*«****'***«*«'*'**'**'**'*****'*«***'***'*'*'***«*******« 

CLOSE  INPUT  AND  OUTPUT  FILES  * 

*********************************************************** 


close (14) 
close (15) 
stop 
end 
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program  convert imp 


lliis  program  %^itten  by  Donald  T.  Hooker  II  on  3/15/93  * 

Revised  3/30/93  * 


*  This  program  converts  numerical  data 
f ran* imp . unconverted*  * 

*  into  ascii  format  for  use  in  modifying  ingrido  for  modal* 

*  imperfections  * 


VARIABLE  DECLARATION 

f*  -  data  to  be  converted  to  ascii  format 
m  ~  nxunber  of  nodes 


character* 13  fl 
character*!  f2,f3,fz 
character*18  f4 
character*!  f5,f6 
character* !8  f7 
character*!  f8,f9 
character*!8  f!0 
character* 5  f!! 

open  (!4,  files '  iitp. unconverted' ) 
open (15/  files  •  inqp.dat ' ) 
fas' O' 


*  READ  NUMBER  OF  NODES  TO  BE  MODIFIED  * 


read (14,*)  m 

*  *************************************************************** 

*  DO  LOOP  TO  READ  DATA  AND  CONVERT  IT  TO  ASCII  FORMAT  * 

*******************************ir*'*lririfitiiii*iiiiifkit*1HHiitii*1t*lt*1t1t* 


do  100  is!,m 

read(!4,998)f!,f2,f3,f4,f5.f6,f7,f8,f9,fl0,f!! 

998  format ( al3 , a! , a! . a!8 , a! , a! , al8 , a! , a! , a!8 , a5 ) 

write(!5, 998)  f !,  f3,  fz,  f4,  f6,  f  z,  f7,  f9,  f  z,  f!0,  f!l 
100  continue 

*********************'*****'******************************i»*** 

*  CLOSE  INPUT  AND  OUTPUT  FILES  * 


close (14) 
close (15) 
stop 
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